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ANALYSIS OF RADIOMETER CALIBRATION EFFECTS WITH TOUCHSTONE 

By 

William D. Stanley* 

ABSTRACT 

The microwave circuit analysis program TOUCHSTONE is used to study two 
effects of importance in radiometer calibration. The two effects are 
impedance mismatches at the antenna- air and cold load- air interfaces and 
dissipatives losses, which radiate thermal noise into the system. The 
results predicted by TOUCHSTONE are shown to be in very close agreement 
with earlier results obtained by purely analytical methods. The techniques 
used in establishing the circuit models and in processing the resulting 
data are described in detail. 


* Chairman, Department of Electrical Engineering Technology, Old Dominion 
University, Norfolk, Virginia 23529-0243. 
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INTRODUCTION 


Two phenomena of considerable interest in radiometer remote- sensing 
systems are those that arise from (a) multiple reflections and (b) 
dissipative losses. Both phenomena may produce errors in either the refer- 
ence calibration process or in an operational mode. Analytical methods are 
very cumbersome to apply and often do not provide an intuitive insight into 
the process involved. 

In an effort to provide an alternate approach to this problem, the 
possibility of computer modeling of the phenomena has been investigated. 

In particular, the program TOUCHSTONE has been used with much success. 
TOUCHSTONE Is a trademark of *EESOF, Inc., and it Is a comprehensive micro - 
wave analysis program providing features for modeling both discrete and 
distributed parameter high frequency and microwave systems. 

Complete operational details of TOUCHSTONE are provided in a user's 
manual , and they are too unwieldy to attempt a summary coverage in this 
report [reference 1]. Instead, the emphasis here will be directed toward 
the development of the circuit forms that may be represented with various 
TOUCHSTONE models and how they may be related to actual radiometer circuit 
behavior. Some data obtained from an actual radiometer system will be 
analyzed, and the results will be compared with an earlier analysis per- 
formed on that particular system. 

Two effects to be delineated In this study are the representation of 
mismatches and the representation of dissipative losses. The first effect 
can be achieved with transformer models, and the second effect can be 

*EESOF, Inc. 

31194 Albania Drive 
Westlike Village, CA 91362 
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achieved with resistors coupled to the system with transformers. The 
various S parameters can be determined and related to the power trans- 
mission coefficients for the signal, the dissipative losses, and system 
mismatches. A radiometer calibration system containing multiple reflec- 
tions will be analyzed in detail. 


SIMULATION OF A MISMATCH 

All measurements with TOUCHSTONE are based on the assumption that each 
port is terminated in a reference resistive characteristic impedance. The 
default impedance is Z Q - R 0 - 50 0, although this value may be changed if 
necessary. 

A desired mismatch may be introduced at any port by using a trans- 
former to couple between the port and the remainder of the system as illus- 
trated in Figure 1. Let T represent the desired reflection coefficient at 
the given port, and assume a reference characteristic impedance of 50 0. 

The input impedance Z^ n as seen at the load must be such that 

50 - Z, 


'in 


( 1 ) 


50 + Z, 
in 

This leads to 


Z. - 50 (1-r) 

in 


( 2 ) 


i + r 

Assuming that the characteristic impedance of any component connected 
to the right-hand side of the transformer is 50 Q, the turns ratio N must 
satisfy 

.2 


50N 


z. - 50 (l-r) 

in 


i + r 


(3) 
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which leads to 



Either the value given by (4) or its reciprocal will be satisfactory. The 
latter case refers to a reversal in the sign of the reflection coefficient, 
but the magnitude is the same. 


SIMULATION OF DISSIPATIVE UOSS 

Dissipative loss in a transmission system may be modeled with a sym- 
metrical T-pad attenuator. Assume under matched conditions that the net- 
work power ratio may be defined as 

?ln » e 2a (5) 

P . 
out 

where P^ n and P 0 ut are t ^ xe respective input and output values and a is the 

attenuation constant. It can be shown with image parameter theory that 

this attenuation can be simulated with the circuit of Figure 2 . 

The available thermal noise power P av of any resistance R is given by 

P - kTB 
av 


where k - Boltzmann's constant - 1.38 x 10*^3 J/K, T is the absolute tem- 
perature in Kelvin, and B is the bandwidth in hertz. When the input and 
output ports of the circuit are terminated in the characteristic or image 


impedance, the output noise power P no is given by 


P - G kTB - G P 
no av av av 


(7) 


where G av is the available gain relating the maximum available output power 
to the power available from the resistance. 
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Next, consider the T-pad attenuator of Figure 3(a), and note that this 
circuit can be realized with three transformers as shown in (b) . The turns 
ratio of the various transformers must satisfy 


50 


R-, 


N, 


( 8 ) 


and 


50 


- R, 


N„ 


(9) 



( 10 ) 


(ID 


Denote the three ports corresponding to the three resistors with the 
subscripts " 1 ", " 2 ", and " 3 ", and denote the input and output ports with 
the subscripts "i" and "o", respectively. Let Sqi, Sq 2 > and s 03> represent 
the scattering parameters relating the transmission between ports 1 , 2 , and 
3, and the output. The net available output noise power P no is given by 


no 


- IS 


oil 


kT l B+ s 02 


2 kT 2 B+|S Q 3 | 2 kT 3 B 


where T^, T 2 , and T 3 are the temperatures of the three resistances, respec- 
tively. The net output noise temperature T no resulting from the three 
resistances can be expressed as 


4 


(13) 


T 

no 


P 

no 


kB 
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Next, assume that T]_ - T 2 “ T 3 - Tp, in which case 


no 


1 1 2 i 2 2 

I S 01 1 + i s 02 I + l s 03 


(14a) 


- dpT 


(14b) 


where 


l d | 2 “ I S 01 I + I S 02 I + I S 03 I 


(15) 


The quantity d 2 is the dissipation factor for the network, and it is the 


weighting factor that determines the contribution to the output noise 
through loss from the physical temperature Tp* From the preceding analy- 
sis, it can be deduced that a knowledge of Sqi* S()2> an ^ ^03 ma y deter - 
mined with TOUCHSTONE. Further, the use of the processing function permits 
the square of the magnitude to be determined within the program. 


REDUCTION IN NUMBER OF ATTENUATOR PORTS 
Realization of the attenuation network with three resistors obviously 
requires three separate ports in addition to input and output ports. A 
procedure will now be developed for eliminating one of the ports by means 
of an impedance transformation of the original circuit. 

Referring to the T-pad of Figure 3(a), the open circuit impedance 
matrix [Z] can be expressed as 
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R 1 + R 2 


[Z] - 


R 1 + R 2 


(16) 


The impedance level of the network at the input can be multiplied by a 
constant 1/N multiplying both the first row and the first column by 1/N . 
The modified matrix [Z] ' is then expressed as 


[Z]' 


(R^ + R2) ^2 


N 


2 


N 



R, + 


(17) 


Assume next that it is desired to realize [Z]' as an L-network with the 
first element on the left as the series element. For this to be possible, 
the following equality must be satisfied: 


N 


R 1 + 


(18) 


or 


N 


R 1 +R 2 


(19) 


To compensate, an ideal transformer with turns ratio N/l from primary to 
secondary must be placed at the input. The resulting circuit, after calcu- 
lation of the two resistance values, is shown in Figure 4. 
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To simulate the circuit with TOUCHSTONE, two additional transformers 
are required to transform the 50-fl terminations to the values of the two 
resistances. The resulting form is shown in Figure 5. The values of N^, 
Ng, and Nq must satisfy 


N. 


R„ 


R 1 + R 2 


50 

N 


B 


R 1 (R 1 + V (R 1 + 2R 2 ) 


R„ 


or 


n b - 


^50 1 


Vr i (R 1 + R 2 ) (Rj_ + 2R 2 ) 


50 


N, 


R 1 + R 2 


( 20 ) 


( 21 ) 


( 22 ) 


(23) 


isj _ ^ (24) 

w c ~ 

>/ R 1 + R 2 

Although many different attenuation values were considered in this 
study, the preceding process is illustrated for a representative case in 
Figure 6. A T-pad attenuator producing an insertion loss of 1 dB is shown 
in (a) . The circuit of (b) represents a TOUCHSTONE model employing three 
dissipation ports. By employing the transformation previously described, 
the circuit is reduced to one having two resistances as shown in (c) . 
Finally, the TOUCHSTONE realization is shown in (d) . 
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CALIBRATION SYSTEM 


The basic system configuration for calibration of a noise- injection 
radiometer is shown in Figure 7. A cold load at temperature is located 
at a distance d from the radiometer. In most cases of interest, this load 
is cooled to the temperature of liquid nitrogen, in which case Tl — 77 K. 

A closed- loop feedback noise-balancing Dicke radiometer has been 
assumed. Injected noise is added to the input noise such that the mean 
noise power levels on the two halves of the Dicke Cycle are equal. The 
physical temperature of the radiometer is denoted as Tp, of which a typical 
value is T p - 308 K. Various transmission line models of the calibration 
system have been considered in an earlier report [reference 2], Equations 
were derived for representing the combined effects of the mismatches at the 
cold load-air interface and the antenna-air interface. In this report, 
these models will be simulated with TOUCHSTONE. 

SIMULATION CASES 

In order to provide a comparison with results obtained through purely 
analytical methods, some actual data considered in an earlier report were 
used [reference 2] . The values assumed are 

Antenna Reflection Coefficient - 0.141 (-17dB) 

Cold Load Reflection Coefficient — 0.05 (-26dB) 

Cold Load Temperature - 77 K 
Reference Physical Temperature - 308 K 
Several different radiometer configurations were considered. The 
different configurations, along with the TOUCHSTONE file names, are as 
follows : 
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RADMET1 : lossless system without a circulator 

RADMET2 : lossless system with a circulator 

RADMET3 : lossy system without a circulator 

RADMET4: lossy system with a circulator 

In an actual radiometer, the variation of measured temperature is 
normally observed as a spatial function of the distance between the load 
and the antenna. With TOUCHSTONE, however, it is easier to maintain a 
constant spacing and vary the frequency since the observed results are the 
same. The path length was established as a one-eighth wavelength at a 
frequency at 1 GHz . 

As a comparison of the different models, various runs were obtained, 
and certain variables were determined and tabulated. The runs are summa- 
rized in Table 1. The following quantities are defined: 

1 1 1 ^ — power transmission coefficient from Input to output 

(or transducer gain) 



power transmission coefficient from input to output 


2 


power dissipation function 



net output power reflection coefficient 


The net output temperature T out is given by 


out 


mVkiVmVw 2 ^ 


(25) 


where T ou ^ is the net output temperature, Tl is the load temperature, is 
the circulator temperature, Tp is the physical temperature of the dissi- 
pation losses and Tp is the reverse temperature radiated back toward the 
antenna junction by the radiometer. 
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The variations In each of the power coefficients were observed over a 
frequency range encompassing at least one full cycle. The maximum and 
minimum values of these quantities are denoted with the subscripts "max" 
and "min", respectively. The first configuration considered was a lossless 
system without a circulator (RADMET1) , and the layout In shown in Figure 8 . 
Since there is neither a circulator nor dissipative loss, no data appear 
under those applicable columns. It turns out that the maximum value of the 
transmission coefficient occurs at the same point as the minimum value of 
the output reflection coefficient and vice-versa. 

The second configuration considered was a lossless system with a 
circulator (RADMET2) , and the layout is shown in Figure 9. The maximum and 
minimum values of the pertinent data are shown in the second row of 
Table 1. 

The presence of the circulator results in y^ ~ 0; i.e., the reverse 

radiation is diverted into the circulator terminating resistor. However, 

noise generated by the circulator termination propagates toward the antenna 

and Is reflected back into the radiometer. Consequently, the values of 

It I 2 and ItJ 2 , , in this case are exactly the same as for 

I Cl max I Cl min 

9 i 1 2 

| -y | , and 7 , . for the preceding case. 

I 'o I max I o I min 

A conclusion that can be drawn from the preceding discussion is that 
there is no advantage gained from placing a circulator in the front end of 
the radiometer if the circulator termination is at the same physical tem- 
perature as the radiometer. Indeed, if there is any loss or mismatch 
produced by the circulator, the system may be further degraded. However, 
if the circulator is nearly ideal, and if its termination can be cooled to 
a lower temperature, some advantage could be gained. 
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The third configuration considered was a lossy system without a circu- 
lator (RADMET3) , and the layout is shown in Figure 10. At this particular 

point in the analysis, a loss of 1 dB is assumed. (Later, the results of 

varying the loss will be considered.) The maximum and minimum values of 
the data obtained are shown in the third row of Table 1. The 

I 1 2 

presence of dissipative losses results in a non-zero value for |d| . It 

2 

turns out that the maximum value of |t| corresponds to minimum values of 

2 i 1 2 

| d | and |y o | and vice versa. 

The fourth configuration considered was a lossy system with a circu- 
lator (RADMET4) , and the layout is shown in Figure 11. As in the preceding 

case, a loss of 1 dB is assumed at this time. The maximum and minimum 

values of the data obtained are shown in the fourth row of Table 1. As in 
the lossless case, the output reflection coefficient is now zero, but the 
circulator termination now contributes to the output, and the weight of 
this contribution is the same as for the output reflection function when no 
circulator is used. 

DISSIPATIVE LOSS EFFECTS 

A study was made to determine the net output temperature as a function 
of dissipative loss. Although either RADMET3 (without circulator) or 
RADMET4 (with circulator) could be used, the former was chosen because it 
has one less port in the simulation. As a prelude to the analysis, values 
of the various parameters required to simulate the transmission as a func- 
tion of the decibel loss were calculated. The results are tabulated in 
Table 2. 

The results of various simulations with RADMET3 are provided in Table 

2 

3. As indicated earlier, the maximum value of |t| occurs at the minimum 
values of |d|^ and vice versa. Further, the only difference in the 
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results that would be obtained with RADMET4 is that the values listed under 
the 1 7 | ^ columns would be transferred to the 1 | columns. 

After the preceding data were obtained, various temperature parameters 
were calculated and tabulated in Table 4. The definitions of these quanti- 
ties are listed as follows: 


- cold load temperature; i.e., the temperature that 
should be measured 


T ma x - maximum temperature as sensed by radiometer 

X m in - minimum temperature as sensed by radiometer 

f_.„ tm mean value of measured temperature 

UlcaS 

T bias “ fixed bias error in measurement 

T var - peak value of variable component of error in measure- 
ment 

Formulas for determining the preceding quantities are the following: 


T + T . 
max min 


(25) 


me as 


T, . - T - 77 

bias meas 


(26) 


T - T , 
max min 


(27) 


var 


A comparison between all points common to this study and the earlier 
analytical study [reference 2] is given in Table 5. The results of the 
earlier study were based on a few simplifying assumptions in the final form 
of the equation, so perfect agreement is not expected. However, it is 
obvious that the results are very close to each other and differ by less 
than 1 % in all cases. 
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TABLES 
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LOSS 

n 

V n 

NA 

NB 

NC 

0 dB 

0 

oo 

i 

oo 

0 

0.1 dB 

0.2878 

4343 

0.9999 

9.320 

0.1073 

0.2 dB 

0.5756 

2171 

0.9997 

6.589 

0.1517 

0.3 dB 

0.8634 

1447 

0.9994 

5.379 

0.1858 

0.4 dB 

1.1511 

1085 

0.9989 

4.659 

0.2146 

0.5 dB 

1.439 

868.1 

0.9983 

4.163 

0.2398 

0.6 dB 

1.726 

723.2 

0.9976 

3.799 

0.2626 

0.7 dB 

2.014 

619.7 

0.9968 

3.515 

0.2836 

0.8 dB 

2.303 

542.1 

0.9958 

3.284 

0.3031 

0.9 dB 

2.588 

481.7 

0.9947 

3.096 

0.3213 

1.0 dB 

2.875 

433.3 

0.9934 

2.934 

0.3386 

1.5 dB 

4.307 

288.1 

0.9853 

2.383 

0.4135 

2.0 dB 

5.731 

215.2 

0.9741 

2.048 

0.4757 

2.5 dB 

7.146 

171.3 

0.9600 

1.814 

0.5293 

3.0 dB 

8.550 

141.9 

0.9432 

1.636 

0.5765 


Table 2. Attenuation network parameters as a 
function of insertion loss. 
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Table 3. Transmission and dissipation factors for lossy system without 
circulator as a function of insertion loss. 






LOSS i&ax i min max min meas bias 
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Table 4 . Temperature variables for lossy system as a function of insertion loss 


mss 


T (i) 

max 


m(2) 

max 


% Diff 


T (l) 

min 


m(2) 

min 


% 


0 dB 

85.43 

85.32 

0.13 

78.91 

78.93 


0.1 dB 

90.49 

90.21 

0.31 

84.13 

84.32 


0.5 dB 

109.63 

109.88 

0.23 

103.83 

103.73 


1.0 dB 

131.20 

131.36 

0.12 

126.03 

126.10 


2.0 dB 

167.57 

167.55 

0.01 

163.46 

163.52 


3.0 dB 

196.45 

196.43 

0.01 

193.19 

192 . 15 



(^Earlier study [reference 2] 

( 2 ) TOUCHSTONE simulation 

Table 5. Comparison of earlier analytical studies and TOUCHSTONE 
simulation studies at coiranon points. 


Diff 


0.03 

0.23 

0.10 

0.06 

0.06 

0.54 
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MM 


2.875 ft 


2.875 ft 



50ft 50ft 



Figure 6. Representative attenuator circuit (1 dB) converted to 
two-port form. 


cold 

load 
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Figure 7 


antenna directional 



Basic radiometer system configuration for calibration. 
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